Market price simulator based on analog 
electrical circuit 

Aki-Hiro Sato^*and Hideki Takayasu^ 
Department of Applied Mathematics and Physics, 
Kyoto University, Kyoto 606-8501, Japan, and 
^ Sony Computer Science Lab., Takanawa Muse Bldg., 3-14-13, 
Higashi-Gotanda, Shinagawa-ku, Tokyo 141-0022, Japan. 

February 1, 2008 



Abstract 

We constructed an analog electrical circuit which generates fluctuations in 
which probability density function has power law tails. In the circuit fluc- 
tuations with an arbitrary exponent of the power law can be obtained by 
adjusting the resistance. With this low cost circuit the random fluctuations 
which have the similar statistics to foreign exchange rates can be generated 
as fast as an expensive digital computer. 

Key words, random multiplicative process, power law, foreign currency 
rate, analog electrical circuit 

1 Introduction 

Random multiplicative process (RMP) is attracting researchers as a new 
mechanism of generating fat tails in the distribution of price changes in open 
markets recently (Takayasu et al. 1997, Jogi et al. 1998, Sato et al. 2000). 
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One of the mechanisms of the fat tails is a feedback mechanism of both 
positive and negative in a random manner. The positive feedback plays a 
role of amplification while the negative one plays a role of damping. In the 
presence of positive feedback one may consider that the system is unstable, 
however, this is not always correct. When the damping effects dominates the 
amplification the system can generate statistically stationary fluctuations. 

Recently fat tails in the distribution have been confirmed in the market 
price fluctuation (Mantegna et al. 1995). Among several approaches to this 
phenomena (Bak et al. 1997, Lux et al. 1999), the authors demonstrated that 
RMP can be derived from a microscopic artificial dealer's market simulation 
model (Sato et al. 1998). The RMP has also been derived from a macroscopic 
theoretical analysis of the critical dynamics of the balance of demand and 
supply (Takayasu et al. 1999). The multiplicative effect of the market comes 
from the dealer's forecasting which follows the trend of the latest market 
price changes, which obviously enhances the fluctuations in price. 

In this article we review the analog electrical circuit which generates 
power law fluctuations based on the RMP theory. The application of this 
circuit is to generate fluctuation which is statistically similar to market price 
change in inexpensive way at high speed. The goal of this study is a risk 
estimator that calculates realizations of market prices directly by parallel 
processing units and calculates a statistical feature of the fluctuation from 
these simulated data in the meaning of an ensemble (see. fig. |l].). 

The outline of this article is the followings. In sec. |] we analyze tick data 
of yen/dollar exchange rates, and show that a probability density function of 
the market price changes has fat tails of roughly power law. In sec. |] a ran- 
dom multiplicative process is formalized by a stochastic differential equation 
with a multiplicative noise. It is indicated that the probability density of a 
dynamical variable follows a power law distribution. In sec. ^ an example 
of electrical circuit diagram is shown and an equivalent equation of the pro- 
posed circuit is described. In sec. ^ we show the result of observation and 
we will discuss them. Sec. |] is devoted to the concluding remarks. 
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2 Statistical properties of foreign currency 
rate 



Mantegna et al. investigated time series of stock market index of S&P500 
and reported in their famous paper (Mantegna et al. 1995) that a probabihty 
density function of changes obeys a power law distribution. They indicated 
that the power law exponent is estimated as 1.4. However, it is clarified that 
the exponent is not universal by other researches. Besides stock market prices 
it is known that the probability density function of foreign exchange rates 
also follows a power law distribution (Takayasu et al. 2000). We analyze 
tick data of yen/dollar exchange rate of 578,508 data points from January 
to March 1999. Let us denote the rate of the sth tick as Vg. The change 
of the rate is defined as Avg = — Vg-i. Fig. |^ shows typical examples 
of time series of yen-dollar exchange rates and corresponding changes. The 
cumulative distribution function (CDF), which is defined by 

/ — \x\ roc 
p{x')dx' + / p{x')dx', (1) 
-oo J\x\ 

is shown in fig. |^. It is clear that the CDF has a liner slope with a quick 
decay by 1 yen in log-log plots. From the hner slope in log-log plots we can 
estimated the exponent of the power law distribution, 

F{> \x\) oc \x\-^. (2) 

The best fit value of f3 we obtain from the exchange rate is (3 = 1.8. 



3 Theory 

In this section we show a brief outline of a continuous time version of ran- 
dom multiplicative process. The random multiplicative process is given by a 
stochastic differential equation, 

-j^ = u{t)v{t)+at), (3) 

where v{t) is a dynamic variable and z/(t) and C,{t) represent a multiplicative 
noise and an additive noise, respectively. Denoting an ensemble average by 
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(...), we assume the following relations for both multiplicative and additive 
noises based on the Gaussian white noise theory. 



Ht)) = u, (4) 

- iy][u{t2) - u]) = 2DJ{h~t2), (5) 

m) = 0, (6) 

iatinh)) = 2D^6{h-t2), (7) 

where u represents an average of the multiplicative noise, and and 
represent the strength of the multiplicative noise and that of the additive 
noise, respectively. 

The probability density function of vit) in eq. (|^) is denoted as p{v,t), 
is known to follow the generalized Fokker-Planck equation (Deutsch 1994, 
Venkataramani et al. 1996, Nakao 1998). 

-p{v, t) = —[-{u + D,)vp{v, t) + Q^{Dy + Di:)p{v, t)] . (8) 

By solving eq. (||) for a steady state §ip{v, t) = and boundary condition, 
^p(±f„, t) = 0, we get the stationary distribution. 

p{v) (x{D^ + Dy)^-^, (9) 

which has power law tails for large f , 

p{v) oc \v\-^-^, (10) 

where 3 = In other words the power law exponent is given by a simple 

function of the average and variance of the multiplicative random noise. 



4 The Circuit 

From the assumption for the multiplicative noise in the RMP theory de- 
scribed in sec. ^ it is important that i'{t) takes both positive and negative 
values to realize the power law tails. With an electrical analog circuit this 
means that it is necessary for the circuit to include both positive and negative 
feedbacks in the meaning of probability. We solve this problem by using an 
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analog multiplier. Our block diagram of circuit and an implemented noise 
generator are shown in fig. ^ Vo in the figure represents an output volt- 
age, and is an output directly from the noise generator. As it is seen 
from this figure the circuit contains the noise generator, an analog multiplier 
(Analog devices, lOMHz, 4-quadrant) and an operational amplifier (National 
Semiconductor, LF157) for integrator (Sato et al. 2000). 

The output of the noise generator plays a role of the multiplicative noise 
in eq. (^. In the noise generator a shot noise of the zener diode in fig. ^ is 
amplified by an operational amplifier, and it is output from generator through 
high pass filter. LF157 has the 20M product of a voltage gain (G) and a 
bandwidth (B). The bandwidth is given by B=200kHz in the noise generator 
because we put G=100. Thus we expect a frequency characteristics of the 
noise generator to be up to 200kHz . Let us explain how to realize both 
positive and negative feedbacks in the meaning of probability. We realize 
both positive and negative feedbacks by multiplying Vo with the output of the 
noise generator /i(t) in the analog multiplier and by connecting the product 
to the negative input of the operational amplifier for integrator. Although 
the additive noise term, ^{t), is not explicitly added in the circuit, it derives 
from either thermal noises of the operational amplifier or from an external 
electro-magnetic noises. 

An equation equivalent to the circuit diagram of fig. ^ is given as 

' ^ ■ ^Mt))vo+m, (11) 



dt ^RfC R^C 

where /c is a factor of the multiplier and k = 1/10. ^(t) represents the 
additive noise effect. The strength of multiplicative noise fi{t) depends on 
the value of a variable resistor in fig. ^ because R^ is a factor of fi{t) in 
eq. (0). Therefore, we expect that the power law exponent for the output 
is a function of Ry. 



5 Results of observation and discussion 

We measure the output of the circuit Vo for R^ since we expect that the power 
law exponent /3 is a function of R^ from the above discussion. We obtained 
the output through a 12-bit AD converter (Microscience, ADM-652AT) and 
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processed it as digital data in digital computer. A sampling frequency is 
125kHz throughout all the observations. 

We show a typical example of time series of Vq at i?„ = 50fi and lOOfi in 
fig. ^. The time series at = 50Q sometimes exhibits larger fluctuations 
than = lOOn. We show log-log plots of CDF oiv{t) at R^ = 25,50,75,1001) 
in fig. ||. Log-log plots of the CDF have liner parts for about a decade. We 
clearly find that the exponent (3 depends on Ry as expected qualitatively. 
Moreover as shown in fig. |^ we find a linear relation between /3 and R^. 
The autocorrelation function R{t) and the volatility autocorrelation function 
R^'^\t) are defined by 



We show the autocorrelation function and volatility autocorrelation function 
in fig. H. Both autocorrelation function and volatility autocorrelation func- 
tion have quick decay. It is known that the volatility autocorrelation function 
of real market price changes has a long time correlation. The disagreement 
of the volatility autocorrelation occurs because the model equation is too 
simple to describe time structure of fluctuations. We need to study higher 
order differential equations with multiplicative noises. 

6 Conclusions 

We proposed an analog electrical circuit as an analog generator of power law 
fluctuations. We described a theoretical equation representing the electrical 
circuit and showed that the probability density function of a dynamical vari- 
able has power law tails. We measured the output of the circuit and observed 
its cumulative distribution functions. The cumulative distribution of output 
voltage has power law tails. The power law exponent can be tuned by con- 
trolling the variable resistance Ry. We expect that the proposed circuit is 
applicable to generate fluctuations having power law distribution in a much 
cheaper way than any digital computing methods. Moreover, fluctuations 
of the circuit may be of use for risk estimation in foreign exchange or stock 
market in the near future. 
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Figure 1: Conceptual illustration of risk estimator that calculates realizations 
of market prices directly by parallel processing units. 
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Figure 2: Typical examples of yen/dollar exchange rates (the above) and its 
changes (the bellow) from January 1999 to March 1999. 
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Figure 3: The log-log plots of the cumulative distribution function of 
yen/dollar exchange rates. A sohd hne represents the power law with P — 1.8. 




Figure 4: Block diagram of the circuit with a noise generator. The circuit 
contains a noise generator, a analog multiplier and an operational amplifier 
for integration. Rf — lOOkQ, C = 10pF,i?.y = 200Q. A variable resistor 
under the operational amplifier is for adjustment of the offset. 
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Figure 5: Typical examples of time series at — 50Q and lOOQ. 
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Figure 6: Log-log plots of CDF of the output Vo{t) at = 25, 50, 75, lOOO. 
Each CDF shows straight lines with different slopes between a decade. 
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Figure 7: The power law exponent l3 plotted against values of a variable 
resistance R^. The power law exponent (3 is numerically estimated from a 
slope of a liner part in log-log plots of CDF. 
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Figure 8: Autocorrelation function and volatility autocorrelation function of 
output Vo{t). 
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